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We use operator product expansion (OPE) techniques to study the spectral functions of currents 
and stress tensors at finite temperature, in the high-energy time-like region uj ^ T. The leading 
corrections to these spectral functions are proportional to ^ T'' expectation values in general, and 
the leading corrections ~ g'^T'^ are calculated at weak coupling, up to an undetermined coefficient 
in the shear viscosity channel. Spectral functions are shown to be infrared safe, in the deeply virtual 
regime, up to order g^T^ . The convergence of (vacuum subtracted) sum rules in the shear and bulk 
viscosity channels is established in QCD to all orders in perturbation theory, though numerically 
significant tails ~ r*/(loga;)'^ are shown to exist in the bulk viscosity channel. We argue that the 
spectral functions of currents and stress tensors in infinitely coupled A/" = 4 super Yang-Mills do not 
receive any medium-dependent power correction. 



PACS numbers: ll.lO.Wx, ll.lO.Jj, 12.38.Mh 



I. INTRODUCTION 

This paper is devoted to applying operator product 
expansion (OPE) [ij techniques to study asymptotics of 
real-time spectral functions p{lo, q) at finite temperature. 
We will obtain the leading corrections in an expansion in 
T /lj at large frequencies lo, with T a characteristic energy 
scale of the medium, to the spectral functions of currents 
{J^), scalar operators (mipip) and stress tensors (in the 
shear and bulk viscosity channels). 

The zero frequency, zero momentum limits of spectral 
functions are related in general to hydrodynamical trans- 
port coefficients, which, for the quark-gluon plasma as 
probed by the Relativistic Heavy Ion Collider 0, have 
been the focus of much recent work (for a more ample dis- 
cussion we refer the reader to [2] and references therein). 
At finite time-like momenta, the vector channel spectral 
function is related to the production rate of lepton pairs 
by the plasma ^3] . 

The OPE techniques employed in this paper, on the 
other hand, give information on the deeply virtual time- 
like region, lu ^ T,\q\. Spectral hmctions in this region 
probe the surrounding medium only on short time and 
length scales, which is why they are related, by the OPE, 
to the expectation values of local operators. The OPE 
analysis thus does not cover the region uj ^ T where the 
thermal corrections to the shape of p are the most im- 
portant and where most of the, e.g., lepton pair emission, 
occurs. 

This paper will cover the following applications. First 
of all, the compact expressions obtained in the OPE 
regime can be used as simple consistency checks on more 
complete calculations. For example, we will prove that 
the leading thermal corrections to spectral functions in 
QCD are proportional in^eneral to T"*, as was observed 
in the early calculations '3^ , and we will give their coeffi- 
cients in certain instances. Also, infrared divergences can 
be fully characterized in the OPE regime: we will show 
that, beginning at order g^T*, but not before, certain 



spectral functions cease to be computable perturbatively 
due to the so-called Linde problem. 

An interesting sum rule was recently proposed by 
Kharzeev and Tuchin and by Karsch, Kharzeev, and 
Tuchin , and used to estimate the QCD bulk viscosity 
near the deconfinement phase transition. We will show 
that their sum rule is sensitive to a numerically impor- 
tant ultraviolet tail, ignored in 5j and Ji], which could 
affect their analysis away from very close to the phase 
transition. 

In the shear viscosity channel we will demonstrate the 
convergence of sum rules in asymptotically free theories. 
Furthermore, we will argue that the left-hand side of such 
sum rules is saturated by a one-loop calculation, which 
however will not be performed here. We will also ob- 
serve the possibility that such sum rules could possess, 
in certain theories, strong ultraviolet tails making them 
discontinuous in the free theory limit {g^ 0). 

Finally we will briefly study spectral functions in the 
strong coupling limit of A/" = 4 super Yang-Mills. We 
will argue that power corrections to spectral functions at 
large virtuality are restricted to polynomial terms in the 
momenta and forbidden for currents and stress tensors, 
generalizing observations of Teaney [3] • 

This paper is organized as follows. After describing our 
formalism in section [TTl we apply it in the weak coupling 
regime to the spectral functions enumerated above, in 
section Hill The physical implications of these results are 
then discussed in section ITVl Finally, in an Appendix [XI 
we reproduce a two-loop diagrammatic calculation which 
confirms the OPE prediction in the vector channel. 



II. FORMALISM 

The physical basis of the (Euchdean) OPE f]\ is the 
separation of scales between that a short-sized probe Aa; 
and that of a typical wavelength ^ T^^ in a medium, 
leading to useful asymptotic expansions in (TAx). In 
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Fourier space for the expectation value of a two-point 
function of operators Oi^2 this gives rise to asymptotic 
expansions at high momenta, suppressing arguments 
other than the frequency: 

GEi2(a;E)~E('^»)^- (2-1) 

n tj 

The locahty of the operators On (which we will assume 
to be Hermitian, without loss of generality) follows from 
T <^ loe- The powers dn are determined by renormaliza- 
tion group equations (RGE), to be reviewed shortly. 

Naively taking twice the imaginary part of the ana- 
lytic continuation of Eq. (|2.1[) to Minkowski frequencies 
uj would yield asymptotics for the spectral function, 

Pi2(^)^E(^«)2I-((Z^)- (2.2) 

with p{lu) the commutator function {[O{lu),O{x—0)]). 
Eq. (|2.2p is "naive" in that the asymptotics of an an- 
alytic function in one complex direction do not, in gen- 
eral, determine its asymptotics along other directions. 
It is nevertheless correct whenever p{u') does admit an 
asymptotic expansion in inverse powers of w (or loga- 
rithms), as will be proved shortly. 

Eq. (|2.2p is the main result of this section. We read it 
as a Minkowski-space version of the OPE: heuristically, 
the locality of its operators is a consequence of the small 
time during which a pair of high-energy particles, cre- 
ated by a high-frequency operator, can travel before it 
is reabsorbed by its complex conjugate. Due to this pic- 
ture, once its coefficients are determined we expect it to 
remain valid even in out-of-equilibrium situations, where 
an Euclidean formulation is not available. 



A. Dispersion relations 

Here we justify the passage from Eq. (|2.ip to Eq. (|2.2p . 
assuming that p{w) admits an asymptotic expansion in 
inverse powers of to and logarithms. The argument is 
based on the dispersive representation of the Euclidean 
correlator^, 

/oo If 
^ , . M ^')- (2.3) 

Note that p is always real (or a Hermitian matrix), 
p = 2ImG. P„(we) is a polynomial in we that is not 
determined by the spectral density. In general Eq. (|2.3p 



At finite temperature, this gives a distinguished analytic con- 
tinuation of Ge from the discrete set of Matsubara frequencies 
at which it is strictly defined ,81 . Substituting o^e = —ii^ with 
u in the upper-half plane, Eq. I|2.3| l always coincides with the 
retarded function G-ii{ui). 



is ultraviolet divergent and a subtracted integral must be 
used, but this does not interfere with the present argu- 
ment. 

The basic point is that, if an upper cut-off < A 
were imposed on the frequency integration in Eq. (j2.3p . 
the resulting Ge('^) would admit an expansion in purely 
integral powers of at large lo. Specifically, for p{lo) 
bounded by |a;|~'''^* (e > 0) at large lo, the fc'^ derivative 
of Eq. (|2.3p with respect to \/lo is shown to vanish at 
LJ — oo. This shows that the asymptotic expansion of p 
directly translates into one for Ge modulo integral terms 
(e.g., terms that are killed by taking derivatives). 

Thus it suffices to match the asymptotic expansions 
of Eq. (|2.ip and Eq. (|2.2p termwise. Two cases must be 
distinguished: non-integral power terms (or powers times 
logarithms), and purely integral powers. 

The former case is dealt with with the dispersive trans- 
form of a power law tail p{uj) oc 

/■°° dcj' 1 1 1 >^ Dnp"-" 

J^^ 27r(a;'— icje) 2 sinTra (— iwe)" lo"- ' 

(2.4) 

where p^ is some infrared cut-off and the sum is an 
asymptotic series with n integers. The coefficients _D„ 
in this sum depend on the details of infrared data (here 
on ojq) but the non- analytic term is a clean reflection of 
the large oj behavior of p, as expected. Its imaginary part 
at real oj produces the right asymptotics. In particular, 
the asymptotics at positive and negative uj can be recon- 
structed independently since their contributions are out 
of phase at real we • Taking derivatives of Eq. (|2.4p with 
respect to a gives identities for integrals with logarithms, 
for which the same argument applies. 

The case of purely integral powers in p, dn integer, is 
special because the Euclidean non-analyticity is a single 
logarithm: it could cancel out between the contributions 
of positive and negative p. Such cancellations would oc- 
cur when dn is even and its contribution to p{ijj) is even 
in w, or when dn is odd and its contribution is also odd. 
In these cases, the Euclidean function would be propor- 
tional to l/(— zwe)" with an imaginary coefficient, e.g. it 
has a "wrong" phase. On the other hand, integral terms 
coming from small frequencies in Eq. (|2.3p . are propor- 
tional to l/(— zwe)" with a real coefficient. Thus the two 
are cleanly separated by their phases, and we conclude 
that in all cases the asymptotics of p can be recovered 
from those of Ge. 

In perturbation theory we do not expect such "wrong 
phase" contributions to the Euclidean OPE, and in any 
case certainly none appears at the relatively low orders to 
which we will be working in this paper. Our corrections 
to spectral functions will come solely from non-analytic 
terms in Ge. 

We now comment the assumption that p(w) admits an 
expansion in inverse powers of w, which we have assumed 
in proving Eq. (|2.2p . Possible violations of it at the non- 
perturbative level (due, for instance, to oscillating terms) 
are discussed in [lO]. However, in perturbation theory we 
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find it hard to see how it could fail, for instance because 
there is no scale to provide an oscillation rate. Thus this 
assumption will be made throughout this paper. 



log(p^/w^). Eq. (|2.8p shows that terms which are poly- 
nomial in p can also contain logarithms of fJ-jjg."^ 



C. RGE in Minkowski signature 



B. Renormalization group equations 

The OPE is based on a systematic separation of in- 
frared and ultraviolet contributions (an interesting dis- 
cussion may be found in [9]). A factorization scale n (for 
us, T < /i ^ w) is introduced, and all vacuum fluctua- 
tions from above this scale are integrated over, leaving 
more infrared and state-dependent fluctuations to be ac- 
counted for by the expectation values of operators. Since 
/i <C w these operators can be taken to be local, in a sys- 
tematic gradient expansion. The restriction to vacuum 
fluctuations ensures that this yields operator relations, 
that is, the OPE holds in any quantum state. 

Omitting Lorentz and internal indices, this yields ex- 
pansions of the form. 



(2.5) 

with TO and g standing for various intrinsic mass scales 
and couplings of the theory and fijj[^ its renormalization 

scale. The renormalized operators O^-'^'' obey the RGE, 
with 7 a matrix of anomalous dimensions, 

0=[M9M + 7]oi''\ (2.6) 



Since we interpret Eq. (|2.2p as a Minkowski-space ver- 
sion of the OPE, it is worth specifying how we mean 
the RGE of operators directly in Minkowski space: we 
mean it to be exactly what it is in the Euclidean OPE, 
namely vacuum fluctuations should be integrated over 
but not state-dependent ones. Provided equivalent regu- 
lators are used, this will reproduce the standard running 
of the operators in Euclidean space. 

Variations on this procedure can easily lead to difficul- 
ties due to transport phenomena (e.g., nonlocal phenom- 
ena), as may be illustrated by a concrete example: that 
of computing the contribution of a fluctuation at scale 
^ gT (with g = ^/AnasT) to the expectation value of a 
local operator in a weakly coupled quark-gluon plasma. 
In this case, it would be natural to fully integrate out the 
scale T, which produces, at the one loop level, the Hard 
Thermal Loop effective theory It is nonlocal with 

support on the classical trajectories of plasma particles. 
The conclusion is that, employing any reasonable proce- 
dure, fully integrating out the scales gT < /i < T in real 
time must convert local operators at the scale T to non- 
local operators at the scale gT. This does not happen, 
however, when only vacuum fluctuations are integrated 
over^, as was assumed in the preceding subsection. 



from which we deduce, in the case that 0[^2 ^^'^ inde- 
pendent of fi (as for currents, which we will exclusively 
study in this paper), the RGE for the OPE coefficients: 



(2.7) 



Assuming the absence of microscopic scales between 
the factorization scale fx and p, the coefficient functions 
can depend only on three scales: the momentum p, the 
factorization scale /i and the renormalization scale fJ-jj[^ 
of the theory. The fj,-^ dependence is determined by a 
RGE 



^Mm« 9g 



Oi(p)02-^i?l2(p)0.- 
(2.8) 

The important point for us will be that the coefficients 
D\2 can only be polynomials in the momenta p. This 
is because correlators at unequal positions are RGE- 
invariant and Fourier transforms can always be made in- 
sensitive to the coincidence limit by taking sufficiently 
many derivatives with respect to momenta. 

Were there no right-hand side to Eq. (|2.8p we would 
conclude, on dimensional grounds, that all logarithms 
in the coefficient functions have to be of the form 



III. ASYMPTOTICS OF SPECTRAL 
FUNCTIONS 

At weak coupling, the OPE coefficients for spectral 
functions, Eq. (|2.2p . are products of Euclidean OPE co- 
efficients and anomalous dimensions, which we now com- 
pute in turn. 



A. Conventions 

For concreteness and simplicity we study the Euclidean 
Yang-Mills theory coupled to np Dirac fermions of the 



^ Examples of such terms are the logarithms in the Green's func- 
tions of free theories. Since they originate from ultraviolet diver- 
gences they do not depend on fi, and their non-analytic behaviors 
~ \og{p'^ / fijjg'^) produce the free theory power tails in the spec- 
tral functions. 

^ Examples of equivalent regulators include a sharp momentum 
cut-off on p and dimensional regularization. 

* A heuristic way to understand this fact is by analogy to the 
situation in the band theory of metals, in which completely filled 
bands do not conduct but only partially filled bands do: the 
field-theoretic vacuum is akin to a completely filled band. 
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FIG. 1: Feynman diagrams giving the leading order OPE 
coefficients for currents and stress tensors. 



same mass m, Se 



+ ^rn)i,,, with 

the covariant momentum. Some important operators (in 
Eudidcan notation) will be 
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(3.1a) 



— Wi Wi ~ [irace partj 



Orr,. = 



(3.1b) 
(3.1c) 



T^" + T^" is the traceless part of the full stress tensor. 
We will also consider the trace of the stress tensor, 



2 i [fermion terms] , 



(3.2) 



327r2 

where 6o = ("y^A— |?t-f2f) is the leading coefficient of 
the /3-function {(3{as) 



B. Euclidean OPE coefHcients 

At the leading order in perturbation theory the OPE 
of currents — ipi'^^'ipi is given by the first dia- 
gram of Fig. [1] A pedagogical introduction to this sort 
of calculation can be found in (l^ . Up to dimension four 
operators we find: 



1 



-1 -7 



Or. 



+ ^Tf [6^'°'5'''^q^-5'''^q^'q°'-5^'^q''q°'+6^"'q°'q'^] , 

(3.3) 

with 75 = 71727374 and the covariant derivative acting to 
its right on ?/;. We will drop the first term since it does not 
contribute to spectral functions, being scale- invariant. 

In isotropic media we decompose the current correla- 
tor into transverse and longitudinal components: we let 



(g4,0,0,|g|), set 



■(J3j3) 



and employ = —^S'^^Tf'^ for the traceless operators. 



Eq. (|3.3p then gives: 



G^(9) - (3.4a) 



3q2 

Similarly, the OPE of Om (scalar channel) is, 

nQ/3 , 



G*(g)^ 



3m2 4m2 



2 . 4 qaqpTr 

q^ 



(3.4b) 



(3.4c) 



that of the trace anomaly (bulk viscosity channel) is, 

,22 ( A If^Ql^rr^aiy , G^ 



(3.4d) 



and that of the shear mode T^^ of the stress tensor, as- 
suming isotropy, is: 

^ ^ „2 -'2\ rr^ii 1 G^ 



G''(g)-^(<zi-rjT-+g^. (3.4e) 



The fermion contributions to Eqs. (j3.4d[) and (|3.4ep be- 
gin at dimension- 6 and have been dropped. However, we 
must be kept in mind that the OPE, like any two-point 
function, is really defined only modulo contact terms 
(terms purely polynomial in momenta). 



Contact terms 



The OPE coefficients of G^ in Eqs. ((Oel) and ((3^ 
are purely polynomial in q. According to the discussion 
at the end of section III Bl this means we have to decide 
whether the operators get evaluated at the scale /^j^g or 
w; alternatively, contact terms depending only on 
could be freely shifted in and out of the OPE as just 
mentioned. 

In the shear channel Eq. (|3.4ep it turns out that, had 
we computed the full OPE coefficient for general T^^'^T"-'^, 
we would have found that its coefficient is purely polyno- 
mial and non-transverse (e.g., leading to q^T^^'^T"-'^ ^ 0). 
Evaluating this operator at the scale lu would lead to 
a non-transverse spectral function, which is impossible. 
Therefore, the G^ term in Eq. (|3.4ep must be a pure 
contact term that runs with the scale /ij^ and does not 
contribute to spectral functions. This issue is discussed 
in 0. 

A similar ambiguity makes it possible to shift the co- 
efficients of Tg'^ and Tj^ in Eq. (|3.4ep by p-independent 
constants. Although we believe that this could, in prin- 
ciple, be settled by studying the Ward identities obeyed 
by the fuU OPE of T'^'T"''^, as in the above paragraph, 
this will not be done in this paper. This will translate in 
an indeterminacy for our shear channel spectral function. 

The Ward identities are much harder to exploit in the 
bulk channel because the trace T^^ is subleading at weak 
coupling. Thus it seems hard to determine, without an 
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explicit calculation of running coupling effects, at which 
scale g^G^ is to be evaluated in this channel. Such a cal- 
culation will not be attempted in this work, so we will 
only be able to determine the asymptotics in this chan- 
nel modulo g^G^ . The term proportional to q^q^T^j'''' / 
in Eq. (I3.4dp is unambigous, however, since the Lorentz 
covariance of the OPE forbids the addition of spin-2 con- 
tact terms to spin-0 operator products. 



find the leading (dimension-four) thermal corrections: 

16q!s ql+'q' 



Sp'iq) 



r 



9<?2 

3q^ q^ 
'^asm?G-p 



[2CFrf -npTFT™] , (3.7a) 



[2CFr}'"-nFTFr; 
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2" 



o 



(3.7b) 



(3.7c) 



D. Anomalous dimensions 



5p^{q) 



%ai q^q^ 



2CFrr-(nFrF-i--&o)rr 



The anomalous dimensions matrix of dimension-four, 
spin-two operators reads p^ . acting on the basis 



-G- 



(3.7d) 



Sp'^iq) - ^b-H^ [2C7FrO«-nFTFT™]. (3.7e) 



^=3^ 



2nFrF -4Cf 
-27ifTf 4Cf 



(3.5) 



giving, e.g., the fi dependence of T^'^ in Eq. 



as log 
3^ 



2CfT^"'^ 



riFTFTr 



up to terms of order 



(3.6) 



The combination Om does not run in perturbation the- 
ory, but the bare mass m has ^-dependence given by 
7m = 3asCF/27r -|- ©(ckg), which must be included wher- 
ever it appears explicitly. 

In theories with massless quarks, G^ docs not run at 
one-loop (at this order G^^T'^'- which does not run to 
any order), though explicit powers of as do run according 
to 7a, = -/3(as)/as = 6oas/27r. 



E. Spectral functions 



Frequency-dependent logarithms in Ge are determined 
by the RGE of the operators entering the right-hand side 
of the OPE, Eqs. \og{l/p'^) log(l/w|). To ob- 

tain p we take twice the imaginary part at real lo (or 
employ Eq. ^T^), log(l/w|) 27r. To help clarify 
our conventions we recall the leading order, massless, 

/ X I vac 

p w\ = 

P^{l)\ 

This way we 



zero-temperature results, with q 

J / \ I vac 



47? 



SOtt 



[dA + ^npdF] 0. 



The presently undetermined coefficients C and D are 
due to the ambiguities discussed in subsection IIII C[ with 
G=D=1 corresponding to choosing the renormalization 
scale u! in Eqs. (|3.4dp and (|3.4e[) . G and D are in principle 
both computable by more accurate calculations. For sim- 
plicity, Eq. (j3.7dp accounts for the running of G^ only for 
massless fermions, in general this operator can mix with 
Om- Eqs. (|3.7c|) and (' (j3.7dp ) do not assume isotropy and 
the anisotropic case of Eqs. (|3.7ap and (|3.7bp may be ob- 
tained starting from Eq. p.3p . We note that 6p''' = 
when n-p = 0. 

The operators in Eq. (|3.7p are all Minkowski-signature 
operators, with the energy density being £ = T*"^ = 

2^44 Please also note that Eqs. (|3.7p arc written in 

(H ) metric convention, so both q^ and T''^ have 

opposite sign relative to there Euclidean cousins (e.g., 
T^^ = £-37^ in Eq. ^l^). In Minkowski space with 

fermionic action J2i'^i(^~ Om = V'lV'i- 

At the Stefan-Boltzmann (free) level. 



(T0°> = 

{Om) = 



15 
60 



(3.8a) 
(3.8b) 
(3.8c) 



Eqs. (|3.7[) . together with (|3.8p . constitute the main re- 
sults of this paper. 



IV. DISCUSSION 

A. The photon spectral function 

One motivation for doing the present calculations was 
to resolve, in a logically independent way, a discrepancy 
in the literature regarding the asymptotics of Sp'^{E). 
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This (in fact the complete spectral function for lo ^ gT 
at zero spatial momentum) has been computed a long 
time ago [1, to order (two- loop order), and it was 
observed that the correction was proportional to g^T'^ /uj^ 
at large energies uj. In contrast, a more recent calculation 
[iBl instead observed a (?^T^ behavior. 

The OPE analysis presented in this paper makes 
it clear that the dominant thermal effects must scale 
like g^T^/w^, since the lowest dimension of a (gauge- 
invariant) local operator with a nontrivial anomalous di- 
mensions in QCD is 4. In particular, a g^T^ asymptotic 
behavior is forbidden by the absence of local dimension- 2 
operators. Thus we can confirm (at least qualitatively) 
the early findings [3, 15]. It is difficult here, however, be- 
cause of the different techniques employed, to comment 
explicitly on the careful calculation of [Iq.^ 

In Appendix [Al we reproduce the OPE result Eq. ()4.1|) 
by means of a more standard diagrammatic calculation 
in real-time perturbation theory. 

For p"^ at g = we find (upon reinstating the well- 
known r = result, not computed here): 



47r 



9 



(4.1) 

Interestingly, the correction, though parametrically small 
~ T'^/to'^ at high energies, has a large numerical prcfac- 
tor, suggesting that it could be useful down to not-so- 
large frequencies w ^ T. Comparison with the complete 
two-loop calculation d, [l^ should allow a precise deter- 
mination of the regime of applicability of the OPE, which 
will not be pursued here. 



C. Infrared divergences 

The OPE analysis completely determines the cancel- 
lation pattern of infrared divergences at high energies: 
the OPE coefficients contain only infrared safe zero- 
temperature physics. Infrared sensitivity can enter only 
through the expectation values of the local operators. 

The leading order corrections to spectral functions p 
are proportional to anomalous dimensions times tree- 
level expectation values of dimension-four operators. The 
perturbative series for such expectation values should be 
similar to that for the thermodynamic pressure ^20] . This 
implies that sensitivity of p to the gT scale (and the ne- 
cessity for Hard Thermal Loop resummation) will first 
enter at order g^T'^/oj'^, and that nonpcrturbative physics 
associated with the g^T scale (the so-called Linde prob- 
lem [2l|) will begin to contribute at order g^T'^/iv'^. 

This situation should be contrasted with that for the 
shape of the spectral function at uj ^ T, for which in- 
frared sensitivity shows up much earlier. At ^ T, gT- 
scale physics appears at order g^ |22| , whereas at ~ gT 
it is important already at the 0(1) level (23]. Various 
integrals over the spectral functions, in the spirit of the 
famous T = sum rules [IJ, however, are still governed 
by the Euclidean OPE, which becomes sensitive to the 
gT scale only at order g^ and to the g^T scale at order 
g^ . This seems to constrain, to a large extent, the cor- 
rections to the shape of p to only move things around in 
frequency space. 



D. Bulk viscosity sum rules 



B. 



Massive fermions (m S> T) 



Spectral functions of massive particles (with m ^ T) 
at g = have been considered in [l^. In this case 
the fermionic condensates T^"' and Om do not con- 
tribute, and our results Eqs. p.7p for the thermal cor- 
rections (e.g., ignoring 0{as) vacuum corrections) be- 



come p^\q) ^ ^ 



1 



45 q' 



1 



45 



and p^{q) ~ 
in complete agreement 



with the findings^ of [K 



A priori, the C(l) power tails in G?, Eq. jHjdl, might 
appear sufficiently strong to make even the difference 
(5Gg = Gg— Gg™'^, between the thermal Euclidean cor- 
relator of T^^ and its vacuum limit, require a subtracted 
dispersive representation. However, for asymptotically 
free theories, the RGE invariance of G^ forces us to eval- 
uate the factors of g^ in Eq. (j3.7dp at the scale w, in which 
case p{ijj) ^1/ log^ UJ at worse and an un-subtracted dis- 
persion relation for (5Ge(0) converges: 



duo' 



2ttoj 



-5p^{oj'). 



(4.2) 



^ We do nevertheless agree with the main conclusion ofJl6| , which 
was to rule out the infrared divergences claimed in [iTJ (see also 
the reply [H). 

® To extract the complete asymptotics of we observe that the 
mass subtraction (4.3) must be undone from their final results 
(4.7) and (C.ll). In their notation, at order this means 

adding '^""^^jr'^^ to (4.7) and '^•'^^j'''\ -T-^ + 2M'^T^ /lo^) 
to (C.ll). 



Higher order corrections to OPE coefRcients will be sup- 
pressed by powers of 5^(w') ^1/ logw' and will not affect 
convergence. We are not making any assumption here 
about the value of the coupling constant at the scale T, 
only the scale lo' is important to the OPE coefRcients. 

On the other hand, in [25] the Euclidean correla- 
tor (5Gg(0) = limq^o liniw^o (^Gr (w, q) is evaluated, by 
means of broken scale invariance Ward identities ("low 
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energy theorems" ) [25 



SGiiO) 



{£ + V) 



Ai£-3V), (4.3) 



with £ = T°" and P = T". 

Eq. (|4.3p together with the exact sum rule Eq. 
were used recently in 5] and [6|] to obtain information on 
the bulk viscosity C = limi^^o p{uj)/uj near the QCD 
phase transition. 

It is not our goal here to discuss the equality of the left- 
hand sides of Eqs. (|4.2p and (|4.3p nor possible contact 
terms to be added to the right-hand side of Eq. (14.21) : 
this is discussed further by Romatschke and Son [27j . 
However, since the results of [5j and Q were based on the 
assumption that Eq. (|4.2p is saturated by low uj (together 
with an Ansatz for the shape of the spectral function) 
which is clearly in tension with the existence of the tail 
Eq. (|3.7dp . we would like to investigate the importance 
of this tail. 

Let us thus try to estimate the contribution from the 



ultraviolet region uj 



in the pure glue theory. Set- 



ting as{uj) = 27r/&olog(w/AQCD) in Eq. (|3.7d|) yields: 



uv 



C 



°° duj_ 

'mil, "''-^ 

{W~£) 



-Ct:T\ 



27rr' 



00 



log' 



Aqcd 

£ 



Aqcd 



log 



Aqcd 



log' 



(4.4) 



The logarithms in Eq. (|4.4p are never particularly large: 
setting LJniin = SttT we estimate l/log(ci;niin/AQCD) ~ 
&oas(27rr)/27r « 0.4 with 6o = 9 in np = 3 QCD and 
tts = 0.3. Thus we conclude that, at least in the pure glue 
theory, whenever Eq. (|4.3p is not parametrically large 
(e.g., l/Cg large) compared to 0.4£, the sum rule Eq. (|4.2p 
is very much sensitive to the ultraviolet tail and is not a 
clean probe of the uj < uj^w, region. It seems that this 
could affect the analysis of fs"! and ^ , at least away from 
very close to the phase transition. 

The weak-coupling limit of Eq. (|4.2p is particularly in- 
teresting: its left-hand side is of order g^T"^ [g being 
evaluated at the scale T from now on) whereas its right- 
hand side receives a contribution of order g'^T^ from the 
UJ r^T region [H, 



duj ^ , ^ duj dxblaluj^ ,UJ, dpJilalT'^ 

— 5p{uj) « / ' 3 nTi{-) = ^- . 

-p TTUJ Jq ttcj D47r'^ 2 50 

(4.5) 

At order g^T^ there is another contribution: from the 
~ g^T* ultraviolet tail integrated over a ^ l/i?^ loga- 
rithmic range. In the weak coupling limit this contribu- 
tion is well-separated and is given by Eq. (|4.4p evaluated 
at Wmin T . The dominant term is the second term, 
— £'&QQ;s/47r^, which, remarkably, using Eq. (|3.8ap is seen 
to exactly cancel Eq. (j4.5p . The presently undetermined 



coefficient C would only become important at order g^T^. 
Thus, the sum rule Eq. (|4.2p is obeyed at order g'^T'^ but 
only when the ultraviolet tail is included. This resolves 
a puzzle raised in [26j . 



E. Shear viscosity channel and discontinuity at 



From Eq. (|3.7ep . the asymptotics of the thermal cor- 
rection 5p^ are proportional to g'^ times an operator of 
strictly positive anomalous dimension 7. Upon resum- 
mation of logarithms its behavior will be proportional 
to {\oguj)~^~°' (the 1 coming from the running of g^ in 
an asymptotically free theory and a > coming from 
the anomalous dimension): thus the unsubtracted dis- 
persive integral Eq. (|2.3p converges. As in the preceding 
subsection, higher order corrections will not modify this 
result because the theory is asymptotically free. The in- 
tegrals should also converge in conformal theories such as 
A/" = 4 super Yang-Mills, because in this case all tails are 
associated with operators of strictly positive anomalous 
dimensions 7, and decay like uj . 

Since the dispersive integral vanishes at we — * 00 we 
can write in general: 



5Gl{uj^) = 5Gl{oo) + 



duj' 



-co 2Ti{uj'-iuj-E) 



Sp'^iuj'). (4.6) 



Note that convergence of the integral implies that 6G^ 
approaches a constant at infinity, so that (5Gg(oo) is 
well-defined. According to [131 , the left-hand side at 
uje^O is determined by hydrodynamical considerations. 
In asymptotically free theories we believe that the OPE 
coefficients of SG^{oo) are saturated by a one-loop com- 
putation.'' Knowledge of both of these ingredients should 
yield interesting exact sum rules, involving, at most, the 
pressure, energy density and chiral condensates of QCD. 
We hope to return to this question in the future. 

Here we wish only to discuss a possible discontinuity 
in the free theory limit g^ ^ oi individual terms on the 
right-hand side of Eq. (|4.6p . if the undetermined coeffi- 
cient D in Eq. p.7ep turns out to be nonzero. Consider 
the term 1"°" o n the right-hand side of the Euclidean 
OPE Eq. ()3.4ep . At g = and w = cxo it contributes 
a finite amount to SGe{oo), which discontinuously 
changes to (T^"" + T°") /(I + npTp /2Cf) at any small but 
finite coupling due to running, Eq. (j3.5p . On the other 



Because other twist-two operators acquire positive anomalous 
dimensions, only three operators can appear in 5Gg(oo): the 
traceless and trace part of the full stress tensor T^", and Om- 
The coefficient of T^^/g^ vanishes at tree level but a one- 
loop computation is needed to find that of T'^ ^ . The coefficients 
of twist-two operators and of Om are determined at the tree 
level, but we believe a one-loop anomalous dimension matrix is 
necessary to carefully separate the total T*^" from other twist- 
two operators, and Om from the trace T^"^. 



hand, at any finite but small coupling one has the 0{g'^) 
tail Eq. (|3.7ep in the spectral function, which is to be 
integrated over a 0{l/g'^) logarithmic range similarly to 
the preceding subsection. Its contribution is thus 0(1). 
It is easy to convince oneself that it exactly compensates 
for the discontinuity in Ge(oo). 

Thus it might happen that equations such as Eq. (|4.6p 
are only continuous at g^=0 when both terms on the 
right-hand side are included. We hope to return in future 
work to shear channel sum rules in QCD and in other 
theories, and to the question of whether they actually 
contain strong ultraviolet tails. In pure Yang-Mills, the 
coefficient D is irrelevant and at the leading order such 
tails are absent. 



F. Strongly coupled Af — 4 super Yang-Mills 

The operator spectrum of strongly coupled multicolor 
Af — 4 super Yang-Mills (SYM) has the very peculiar 
property, that the only operators which do not develop 
large anomalous dimensions ^ A^/^ are protected by su- 
persymmetry and have strictly vanishing anomalous di- 
mensions [28j, where A — g^Nc is the 't Hooft coupling. 
There are no small nontrivial anomalous dimensions in 
this theory, even at finite A. 

One way to find power corrections in spectral functions 
is, as discussed in subsection III A\ if the Euclidean OPE 
coefficients have "wrong" phases. This would certainly 
seem peculiar but we have no general argument against 
this possibility. However, the results of [23| suggest that 
the OPE coefficients of protected operators in A/" = 4 
SYM are identical to those of the free theory, for which 
this definitely does not happen. 

Thus we will assume that power corrections to spectral 
functions at high frequencies are associated with non- 
analytic terms in Euclidean correlators. Without anoma- 
lous dimensions the only sort of non-analyticity allowed 
by the RGE (see subsection III B[) are polynomial terms 
in momenta, times single logarithms \og{p'^ / fj.y^'^) . They 
lead to strictly polynomial terms in spectral functions. 

By dimensional analysis and transversality these are 
strictly forbidden in the spectral functions of currents 
and stress tensors (except if they multiply the unit oper- 
ator). Thus these spectral functions are strictly protected 
against medium-dependent power corrections. It would 
be interesting to see whether polynomial corrections ac- 
tually occur in other spectral functions. At finite but 
large A, we expect power tails ^ lu^''^ with n ~ \^/^. 

Thermal corrections to the spectral functions of R- 
currents and stress_tensors at strong coupling have been 
studied by Teaney ■ 7] and observed, remarkably, to decay 
exponentially fast at high energies. This section general- 
izes his observation. 
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FIG. 2: Real-time diagrams of first topology contributing to 
(with the complex conjugate diagrams omitted). The 
arrows show the time flow along retarded propagators, not the 
charge flow; the doubly-dashed propagator is the fluctuation 

function Grr- 
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APPENDIX A: DIAGRAMMATIC EVALUATION 

OF p{; 

This Appendix reproduces a calculation of the trace 
p^{p) of the current spectral function, using real-time 
Feynman diagrams. We work in the high energy limit 
where p = TP" and drop all terms suppressed by Boltz- 
mann factors oc exp(— p''/2T), but keep all power correc- 
tions. The aim is to confirm the OPE result, Eq. (|3.7p . 
for —11^^ = 2p^ + . For notational simplicity in this 
section we assume n-p — ^■ 



1. Outline of calculation 

The two-loop diagrams contributing to the Wightman 
self-energy Ii^{q) are shown in figs.[2]and[31 The use the 
cutting rule of Weldon [30] , which expresses the Wight- 
man self-energy as a product of two retarded amplitudes 
separated by Wightman (on-shell) propagators (depicted 
as the main cut in the figures). Its physical interpreta- 
tion is as follows: the main cut sums over intermediate 
states, as is expected for a Wightman function, and the 
amplitudes are retarded because intermediate states are 
expanded in a basis of "in" states (e.g., free theory states 
defined at i — > — oo). 

To evaluate the retarded amplitudes on each side of 
the cut we use the so-called Schwinger-Keldysh (ra) for- 
malism, as described in fsTl The resulting expressions 



Alternatively, these amplitudes are the analytic continuation 
of the Euclidean ones |32( . Either way their evaluation at n- 
loop involves no more than n statistical factors. This may be 
contrasted with the rules (of Kobes and Semenoff ) employed in 
the second of reference 3 and PI, in which terms in which all 
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(a) 



(b) 



(c) 



FIG. 3: Real-time diagrams of second topology contributing 
to n^, in the notation of Fig. (2] Not shown, the complex 
conjugates to (b)-(c). 



are summarized graphically in the figures: the arrowed 
propagators are retarded propagators and the propaga- 
tors with the double cut are the fluctuation functions 
of this formalism {Grr propagators); vertices are as in 
ordinary zero-temperature perturbation theory (e.g., no 
complex conjugation appears). 

Following the OPE philosophy we look for propagators 
which can become soft, q ^ T p. Visual inspection re- 
veals that no more than one propagator can ever become 
soft simultaneously: at least two hard particles must tra- 
verse the main cut, and to channel their hard momenta 
to the external legs in all cases requires at least two other 
hard propagators. Thus we will organize the calculation 
around the propagator which becomes soft. 

At order as there is no need for HTL resummation and 
the retarded propagators are temperature- independent. 
The temperature dependence is due to the statistical fac- 
tors entering the Wightman and rr propagators, 



Sr{p) 



if 



p'^ -f i(.p^' 



(Ala) 

5G^t':<,„(p) = -<5^''Gb(p), <55<,>,..(p) = -i>GY{p). 

(Alb) 

with the vacuum cuts obeying Gy{p) — 2IieGf!{p)d{p'^). 
We will not use the explicit forms Gb.f{p) = 
27r(5(p^)nB,F(b°|) until the end of the calculation; up 
to then the sole purpose of the Ansatzes Eqs. (jAlbp is 
to simplify polarization sums. Our metric signature is 
(+— )• 



2. Gluon condensate 

First we allow the gluon propagator in diagrams (a) 
and (b) of Fig. [2] to become soft. Upon evaluating the 



Dirac trace this contribution may be written, 

-n>;i(p) 



D 32 



X I 47r^ 



GB(fc) + GB(-fc) 



72/2 

'■k'-pk 



/2/2 




where we have introduced the abbreviations Ji = J J27t)^ ' 
Ik — l—k, Ip — l—p and Ipk = l—p—k, to be used in 
all what follows, and fc <C ^ ~ p is the soft momen- 
tum. It is kinematically impossible for two denomina- 
tors in Eq. (|A2[) to vanish simultaneously and only the 
real part (e.g., principal value) of the propagators con- 
tributes. The ^-integration is Lorentz-covariant and be- 
comes elementary in the rest frame that is singled out by 
the ^-functions. Thus Eq. (jMI) yields: 



GB{k) 



[p^+k 



2\2 



2p-kA 



■In 



p'+k^ 



In 



p2 



1 + E±^ 



/l 


-(^ 






V 


-(^ 


k-A^ 


2 


1+ 








1+ 


—p-k- 


i). 


1 



(A3) 



with A = ^J [p-kY — p^k^- 

To evaluate the diagrams of the second topology. 
Fig. [3l without having to deal with ill-defined expres- 
sions such as 5{P)/P (which would appear in too literal 
an interpretation of diagrams (b)-(c)), we write their sum 
as a discontinuity. 



-n>[;(p) 



D 16 



GB(A:) + GB(-fc) 



27TSill) 



X 2Im 



72/ 7 

I ip-lk 



{P + iePY{ll + iell) 



(A4) 



where we have also included the contribution with the 
self-energy inserted on the lower propagator. In our kine- 
matic regime the poles of the denominators are disjoint 
and occur at positive energies, f' > 0, Z° > 0. The discon- 
tinuity across the squared propagator l/(/^)^ may be con- 
veniently evaluated by integration by parts along Z^9pM , 
yielding the formula^: 



2Im 



(/2+ie/0)2 



^ j5{ll)2.5r) {i+1'p^)f[1), 

(A5) 



for F{1) any function of I regular at 



propagators carry statistical factors appear at intermediate steps 
(only to cancel out in the end). 



^ An alternative way of deriving this result is to treat the self- 
energy insertion as a correction to the external states, in which 
case at this order one gets thermal mass shifts and wave-function 
renormalization factors. See, for instance, llSll. 
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The total imaginary part in Eq. (|A4[) is the sum of that 
from Eq. (|A5p and from that across the l/l'^ propagator, 
2Iml/(/| + iell) = -2n6{ll). Upon performing the I 
integration we obtain 



GB(fc) 



-4 



p-k 



hi 



1- 



(A6) 



and from when the lower fermion propagator becomes 
soft, 



1 



Jk 



k-(p+k) 



{p+kf 



(A9) 



Our final result for the sensitivity to the gluon distri- 
bution in the medium is the sum of Eqs. (jA3p and (jA6[) . 



Fermion condensate 



The total sensitivity of Il ^j^ t o fermions moving in the 
medium is the sum of Eqs. ((XT)) . and (|X9l) and the 
same objects with (fc ^ — fc). 

4. Expansion in 1/p 



Letting the lower fermion propagator become soft in 
Fig. [2] (a) and in its left-right flip, or the rr propagators 
with similar positions in (c) and its conjugate, gives a 
contribution: 



-n>,^(p) 



-64 / GF(fc) / An' 

Jk Jl 

lp-{p+k) l-k 

{p+ky 

2p-k {p^-k^)i2k-p+P) 



ija+kf 



(p+k)^A 



In 



''pk 



1 p-k-A 



(A7) 



where our notation and techniques are as in the previous 
subsection. Contributions in which the upper fermion 
propagators are allowed to become soft are similar, but 
with k replaced by — fc; these will have to be included at 
the end. 

The diagram of Fig. [3] receives a contribution from 
when the upper fermion propagator becomes soft. 



-n>;i(p) 

_ 1 
vr 



2TT5{ll){Plp-k ~ 2lkllp) 



-i2l^jGAk)j^ a^ + ^dO)^(/^+.dg) 



GF(fc) 



hp' 



■ p-k 



p-k+A ' 



A 



In 



-I p-fe-A 
W3 



(A8) 



Each of the contributions Eqs. ((X3|) . (|M)) (AT} 
and (jA9p is free of infrared divergences and is moreover 
local in k. That is, each admits a Taylor expansion in 
positive powers of p-k/p^ and A^, as is readily verified 
from their parity under A — > — A. 

This is the main point of this analysis: divergences and 
non-localities have cancelled out in the sum over cuts, for 
each individual diagram. A Minkowski-space OPE thus 
works for each diagram. Furthermore, upon summing the 
diagrams, we find that the term of order (p^)*^ cancels in 
the l/p^ expansion of the sum of Eqs. (|A3[) and (jA6[) as 
required by gauge invariance, since this would correspond 
to a non-invariant Afj_A^ condensate: gauge invariance 
upon summing diagrams. 

The leading nontrivial term in the expansion arise at 
order 1/p^, 



g'^GpdF 7rp2 



GF{k) 



S{p-kf 



|A^ 



'GB(fc) 



2{p-kf 



p^ 



k* 

oC-) 



(AlO) 



in complete agreement with the OPE result for 
{2p^+p^){p) when np = 1, Eqs. dSJ]) . 
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